Travel time stability in weakly range-dependent sound channels 
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Travel time stability is investigated in environments consisting of a range-independent background sound- 
speed profile on which a highly structured range-dependent perturbation is superimposed. The stability of both 
unconstrained and constrained (eigenray) travel times are considered. Both general theoretical arguments and 
analytical estimates of time spreads suggest that travel time stability is largely controlled by a property to' of 
the background sound speed profile. Here, 2n/(o(I) is the range of a ray double loop and / is the ray action 
variable. Numerical results for both volume scattering by internal waves in deep ocean environments and rough 
surface scattering in upward refracting environments are shown to confirm the expectation that travel time 
stability is largely controlled by co' . 
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I. INTRODUCTION 

Measurements made during the Slice89 propagation 
experiment 1 , performed in the eastern North Pacific Ocean, 
suggested that in that environment the near-axial energy is 
more strongly scattered than the energy corresponding to the 
steeper rays. Similar behavior was subsequently observed in 
measurements made during the AET experiment 2 - 3 , also per- 
formed in the eastern North Pacific Ocean. Motivated in large 
part by these observations, several studies 4,5 ' 6,7 have been car- 
ried out to investigate the dependence of ray stability on en- 
vironmental parameters. All of these studies have focused on 
ray stability in either physical space (depth, range) or phase 
space (depth, angle). The present study extends the earlier 
work by considering the sensitivity of ray travel times to en- 
vironmental parameters. Our focus on travel times is more 
closely linked to the Slice89 and AET observations than the 
earlier sensitivity studies inasmuch as both sets of measure- 
ments were made in depth and time at a fixed range. 

In this study we consider ray motion in environments con- 
sisting of a range-independent background on which a range- 
dependent perturbation, such as that produced by internal 
waves in deep ocean environments, is superimposed. We con- 
sider the influence of the background sound speed structure 
on both unconstrained and constrained (eigenray) measures 
of travel time spreads. Surprisingly, the conclusion of this 
work is that travel time stability is largely controlled by the 
same property 7 of the background sound speed structure that 
controls ray trajectory stability. Although this work was moti- 
vated by the deep ocean measurements mentioned above, the 
results presented apply to a much larger class of problems. To 
illustrate this generality we include numerical simulations of 
ray scattering by a rough surface in upward refracting envi- 
ronments. 

The remainder of the paper is organized as follows. In Sec. 
U]] the equations on which our analysis is based are presented. 



These are the coupled ray/travel time equations written in 
terms of both the usual phase space variables and action-angle 
variables. In SeclTTTltwo simple expressions for unconstrained 
travel time spreads in terms of action-angle variables are de- 
rived (trivially) and shown to be in good agreement with nu- 
merical simulations. In Sec. Qv]an expression for constrained 
(eigenray) time spreads is derived, again using action-angle 
variables, and shown to be in good agreement with numerical 
simulations. The same result was recently derived using a dif- 
ferent argument by Virovlyanskyi4 The combination of the 
results presented in Sees. fTTTI and HVI provide strong evidence 
that, quite generally, travel time spreads are largely controlled 
by the same property co' — dco/dl of the background sound 
speed structure. Here, 2n/co{I) is the range of a ray double 
loop and / is the ray action variable. In Sec. [V]our results 
are summarized and discussed. Two explanations for why CO 1 
controls travel time spreads are given. 



II. BACKGROUND 
A. Theory 

This paper is concerned with the scattering of sound, in the 
geometric limit, by weak inhomogeneities. Consistent with 
these assumptions, our analysis is based on the one-way form 
(cf. e.g. Ref. 8 and references therein) of the ray equations, 



dp dh dz dh 
dr dz' dr dp' 

and travel time equation, 

— = —-h 
dr ^ dr 



where 



h{p,z,r) = -Jc 2 (z,r 



(la) 



(lb) 



(2) 
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Here, r, which is the independent variable, denotes range; z 
is depth; p is vertical ray slowness; T is travel time; and 
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c(z,r) is sound speed. Equations Q constitute a canonical 
Hamiltonian system with h the Hamiltonian, (z,p) the gen- 
eralized coordinate-conjugate momentum pair, and T play- 
ing the role of the mechanical action. It follows from Eqs. 
([^) and dz/dr = tan<f>, where the ray angle <p is measured 
relative to the horizontal, that cp = sin(f>. We shall assume 
the sound speed to be the sum of a range-independent back- 
ground component, C(z), plus a (weak) range-dependent per- 
turbation, 8c(z,r). This allows us to write the Hamiltonian 
as the sum of a range-independent (integrable) component, 

H(PiZ) = — \/C~ 2 {z) ~ P 2 ' pl us a range-dependent (noninte- 
grable) perturbation, Sh(p,z,r) ~ —8c(z,r)/[C 3 (z)H(p,z)]. 

In the background environment, i.e. when 8c = 0, the ray 
motion is naturally described using action-angle variables 
(/,#). The transformed ray equations maintain their Hamil- 
tonian structure with Hamiltonian H(I) = H (p(/,#),z(7, #)), 
namely 



dl 
d7 



dH 



0, 



d# 
~d7 



<?/7 



<3# ' dr dl 
The travel time equation, in turn, reads 



dT 
~d7 



= Ico H 



dr 



(G-M). 



Here, 



G(z,I)=7Cl± d^JC-^) -7/2(7), 



(3a) 



(3b) 



(4) 



where the ± sign applies to ±p, is the (multiply-valued) gen- 
erating function of the canonical transformation that relates 
implicitly the new variables (7, ■&) to the original variables 
(p,z) through 



dG 

dz' 



# = 



dG 
dl ' 



The action variable is defined by 



1 

7= — 
2n 



dzp = — 
n 



dzJC-\z)-H\ 



(5) 



(6) 



estimates presented below is zero. In those cases where this 
contribution is not identically zero, it is very small. 

In range-dependent environments, i.e. with 8c ^ 0, action- 
angle variables can be defined using the same relations as 
in the unperturbed environment 17 . This is in fact possi- 
ble because the range-independent canonical transformation 
(p,z) <—►(/,#) depends solely on the geometry of phase space, 
characterized by the symplectic two-form dp A dz, which is 
the same with 8c = or 8c ^ 0; it does not depend on 
the Hamiltonian itself (cf. e.g. Ref. It is impor- 

tant to realize, however, that in the presence of a range- 
dependent pertubation 8c, the action variable 7 as defined in 
Eq. is not constant anymore because 77 is not a first in- 
tegral of Q. The Hamiltonian in the new variables takes the 
form 77(7) + 8h(I, &,r), where 77(7) =H(p(I, i?),z(7, #)) and 
5/i(/,#,r) = Sh(p(I,-&),z(I,~&),r). The ray and travel time 
equations then read 



d/ 
d~r 



d8h 



d# 

"d7 



= CO 



dSh 

~~dT' 



and 



dr d# _ „- d , 

— =/—-//- <5/i + —G-/#. 

dr dr dr 



(7a) 



(7b) 



In realistic ocean environments Eqs. @ are difficult to 
apply because in such environments the dependence of 8h 
on 7 and ■& is not known explicitly. Because of this, most 
of the results presented below assume that the environment 
is piecewise range-independent. These results follow from 
piece wise application of Eqs. l|3} This procedure eliminates 
difficulties associated with the lack of availability of an ex- 
plicit expression for 8li(I, r). The endpoint correction term 
d(G — W)/dr in (|3j)) will be neglected throughout the re- 
minder of this paper. 

The reason for introducing action-angle variables is that 
these variables provide the most succinct description of the 
underlying ray dynamics. The equations presented in the next 
three sections that make use of action-angle variables can also 
be written in terms of the original phase space variables (p,z), 
but the (p,z) form of the equations are more complicated and 
obscure the important role, described below, played by co' in 
controlling travel time spreads. 



where the loop integral is taken along an isoline of 77 and thus 
z± correspond to the ray upper (+) and lower (— ) turning 
depths. 

Equations (|3^) show that in the background environment 
each ray trajectory forms a closed loop in (/?,z)-space, or 
phase space, on which 7 is a constant. The line integral in 
Eq. (|6ji can also be expressed as an integral over the enclosed 
area in phase space. It follows from Eq. © that 7 = for the 
axial ray and that 7 increases monotonically as axial ray angle 
increases. For rays making multiple loops G advances by 2nl 
and # advances by 2n each time a loop is completed. Fol- 
lowing a ray, the difference G — W makes small oscillations 
about zero. The term d(G — W) / dr gives small endpoint con- 
tributions to T. The fractional contributions from this term to 
T are nonnegligible only at very short range. The contribution 
from this term to several of the analytical travel time spread 



B. Some qualitative features of scattered ray travel time 
distributions 

Some basic features of ray travel times in deep-ocean envi- 
ronments without and with internal-wave-induced scattering 
are shown in Fig. ^ Travel time and ray depth are plotted at 
r = 2 Mm for rays emitted from an axial source in each of the 
two sound channels shown in Fig. [2] The nonscattered ray 
travel times plotted in Fig. ^correspond to rays with launch 
angles <p Q confined to 4° bands. The scattered ray travel times 
plotted in Fig. [^correspond to an ensemble of 200 rays with a 
fixed launch angle, each in the same background sound speed 
structure but with an independent realization of the internal- 
wave-induced perturbation superimposed. (Independent real- 
izations were generated using the same internal wave field by 
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FIG. 1 : (a) Ray travel time and depth at r = 2 Mm in the CO sound 
channel (cf. Fig. 2) without (solid curve) and with (dots) internal- 
wave-induced sound speed perturbations superimposed. The source 
is on the sound channel axis in both cases. In the absence of internal- 
wave-induced perturbations 8° < <p () < 12°, where <p () is ray launch 
angle. In the presence of internal-wave-induced perturbations <p = 
10° in each of 200 realizations of the internal wave field, (b) Same 
as (a) except that the C14 sound speed was used with 10° < <p < 
14° without internal-wave-induced perturbations, and q> = 12° with 
internal- wave-induced perturbations. 



staggering the initial range at which rays were launched with 
Ar = 5 km; the vertical derivative of internal-wave-induced 
sound speed perturbations has a horizontal correlation length 
shorter than 5 km 13 , so this simple procedure ensures statisti- 
cal independence.) Figure^shows two examples of a portion 
of what is commonly referred to as a timefront, consisting of 
many smooth branches that meet at cusps. 

The internal-wave-induced sound speed perturbation used 
to produce Fig. ^ — an d a ll subsequent numerical calculations 
shown in this paper that include internal-wave-induced sound 
speed perturbations — was computed using Eq. (19) of Ref. 
UAl In that expression y and t were set to zero, i.e. a frozen 
vertical slice of an internal wave field was assumed. The 
range-averaged buoyancy frequency profile measured during 
the AET experiment was used. The dimensionless strength 
parameters E and jti were taken to be 6.3 x 10~ 5 and 17.3, re- 
spectively. Horizontal wavenumber and vertical mode number 
cutoffs of 2k km -1 and 30, respectively, were used. The re- 
sulting sound speed perturbation field is highly structured and 
fairly realistically describes a typical deep-ocean midlatitude 
internal-wave-induced sound speed perturbation. 

Figure^shows that internal-wave-induced scattering is pre- 
dominantly along the background timefront. Internal-wave- 
induced scattering also causes a broadening of individual 
branches of the timefront. This broadening occurs even for 
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FIG. 2: Background sound speed profiles used in the numerical work 
presented in Figs. 1, 3, 4, 8, 9 and 10. 

a single realization of the sound speed perturbation field if 
rays with a dense set of launch angles are plotted. In the 
two sections that follow, simple analytic expressions that de- 
scribe time spreads are derived and compared to numerical 
simulations. These expressions describe: (i) the relatively 
large scattering-induced travel time spreads along the time- 
front seen in Fig. ^ as a function of range, without regard 
to the depth or timefront branch on which a scattered ray 
falls; (ii) the relatively large scattering-induced travel time 
spreads of rays whose turning history and final depth are 
fixed, but whose final range is not; and (ii) the relatively 
small scattering-induced broadening of an individual branch 
of the timefront at a fixed depth and range as a function of 
range or ray double loops. We refer to (i) and (ii) as uncon- 
strained travel time spreads and (iii) as a constrained travel 
time spread. The latter is constrained by the imposition of 
an eigenray constraint. Both unconstrained and constrained 
measures of time spreads will be shown below to be largely 
controlled by £»'(/). 

It is not surprising that time spreads in the presence of 
sound speed perturbations are largely controlled by ©' inas- 
much as travel time dispersion in the absence of pertur- 
bations is controlled by o)' . Equation (|3p), with the term 
d(G — W)/dr neglected, can be integrated immediately to 
give 

T(I;r) = [Ico{I)-H{I)]r (8) 
where, from the first of Eqs. ©), / is constant following a 
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FIG. 3: Hamiltonian h vs. range for one of the scattered rays shown 
in Fig. lb in the presence of an internal-wave-induced sound speed 
perturbation field. 



III. UNCONSTRAINED TIME SPREADS 
A. Spreading along the timefront 

In this section we imagine launching rays from a fixed 
source location with a fixed launch angle in an ensemble of 
oceans, each with the same background sound speed struc- 
ture, but with an independent realization of the internal-wave- 
induced perturbation field superimposed. At a fixed range we 
consider the resulting distribution of ray travel time perturba- 
tions without regard to the final ray depth or the ray turning 
point history. Such a distribution is seen in Fig. ^ 

With the assumption that 7 is piecewise constant following 
a ray, it follows from Eqs. Q, neglecting d (G — W) / dr, that 



ray. It follows that 



dT/d7 = Ico'r, (9) 



where the second of Eqs. ®) has been used. This simple ex- 
pression succinctly describes the travel time dispersion seen 
in Fig. ^ in both environments ft)' is negative for all nonax- 
ial rays; because / is a nonnegative monotonically increasing 
function of | <jC> 1 for an axial source, (|9} describes the decrease 
in travel time with increasing |<p | that is shown in Fig.Q 

For one of the rays used to produce Fig. [2 the Hamilto- 
nian h (whose numerical value is identical to that of 77 + 8h, 
although they depend on different variables) as a function of 
range is shown in Fig. [3] It is seen that h vs. r consists 
of a sequence of approximately piecewise-constant segments, 
separated by fairly narrow transition regions. The transition 
regions coincide with each of the rays upper turning points. 
A useful and widely used approximation, the apex approx- 
imation, assumes that the width of each transition region is 
negligibly small. In the action-angle description of ray mo- 
tion that makes use of the apex approximation, 7 is piecewise 
constant following a ray, making jumps A7 of negligibly small 
width A# at each upper turning depth of the ray; between such 
jumps # advances by 2n. A slightly relaxed form of the apex 
approximation in which At? is treated as a small parameter 
will be used in Sec. [TV] Generally, the apex approximation 
works fairly well in typical midlatitude deep ocean environ- 
ments for rays with axial angles greater than about 10°; it is 
usually a poor approximation for rays with axial angles less 
than approximately 5 ° . More generally, the apex approxima- 
tion can be thought of as a special case of a scattering model 
in which 7 is piecewise constant on a sequence of range in- 
tervals of variable extent. Such a model is used in Sees. [TIT] 
andQy] Although the apex approximation is a special case, it 
is useful to focus on this special case when considering con- 
strained spreads because it is often the case 1,2,3 that only the 
steep ray travel time spreads can be measured experimentally. 



AT= fd£/(£)a)'(/(§))A/(£). (10) 
Jo 

Here, AT = T — Tq, where 7b (r) is the travel time of the unper- 
turbed ray whose initial action is 7o; A7 = I — 7o; and 7 and ft)' 
are understood to be piecewise constant functions of r. If the 
product 7<a' sw IqCo'(Iq) over the domain of 7-values assumed 
by the ray between ranges and r, then (1 1 Oft can be approxi- 
mated by 

AT*I Co'(I ) fd§A7(§). (11) 
Jo 

Successive upper turns are typically separated by about 50 km, 
which is large compared to the horizontal correlation length 
of internal-wave-induced sound speed fluctuations, so each 
perturbation to 7 is independent. This leads to ((AI) 2 ) 1 ! 2 ^ 
r l l 2 , where (•) indicates ensemble average, and, in turn, to 
((AT) 2 ) 1 / 2 ~ r 3 / 2 . This r 3 / 2 -dependence was previously de- 
rived by F. Henyey and J. Colosi (personal communication), 
who also found that this dependence is in good agreement 
with simulations. [Our simulations also indicate that when 
{TQ is valid, then ((AT) 2 ) 1 / 2 ~ r 3 / 2 ; our focus, however, is on 
the dependence of JlQi and on the background sound speed 
structure, via 7ft)'. We note, in addition, that if ((A7) 2 } 1 / 2 is 
computed using an ensemble of rays with the same launch an- 
gle then the r-dependence of this quantity is characterized by 
oscillations in r superimposed on the r 3 / 2 trend. The oscilla- 
tions are caused by the cycling between relatively small time 
spreads when distributions are centered near the cusps where 
neighboring timefront branches join and relatively large time 
spreads when distributions are centered near the midpoint of 
timefront branches. At long range when scattered ray distribu- 
tions are very broad, this effect is greatly dimished.] At long 
range the contributions to J10I or ( II Q from ray partial loops 
near the source and receiver are unimportant. Then, because 
2k/ co is the range of a ray double loop, for a ray that has un- 
dergone N apex scattering events and has N complete double 
loops, il It can be further approximated, 

AT* 2.^1(1,-10). (12) 

Equation (I12> provides an explanation for the observation 
(recall Fig. Q that <p = 10° rays in the CO environment are 
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(a) 



(b) 



FIG. 4: Estimates of unconstrained travel time spreads (upper pan- 
els) and the stability parameter a (lower panels) as a function of 
launch angle in background sound channels CO (a) and C14 (b) (cf. 
Fig. 2). To generate the upper plots, ray simulations based on Eqs. 
(1) were performed by launching rays from a source located at the 
sound channel axis with a fixed launch angle in an ensemble of 50 
oceans, each with the same background sound speed structure, but 
with an independent realization of the internal-wave-induced pertur- 
bation field superimposed. 



more strongly scattered than <p = 12° rays in the C14 en- 
vironment; the same internal wave field was used to gener- 
ate the ensemble of scattered rays, so the difference between 
these figures is due to the difference in the background sound 
speed profiles. The <p = 12° rays in the C14 environment 
have a |o'(/o)| -value that is approximately one-fourth that of 
the <p = 10° rays in the CO environment, and this scattering- 
induced time spread is seen to be reduced by a nearly com- 
mensurate amount, consistent with Eq. d!2l >. According to 
that equation, travel time perturbations are the product of the 
amplification factor I to' /to and a term that depends on the his- 
tory of the scattering-induced perturbations to the ray action 
variable. 

It is convenient to introduce the stability parametei&^Jl 



a = Ito'/to, 



(13) 



which, in addition to appearing in Eq. H2\ . is a natural nondi- 
mensional measure of to' '. Figure|4]compares numerically sim- 
ulated time spreads {(AT) 2 ) 1 ! 2 as a function of launch angle 
tp for a source on the sound channel axis at r = 2 Mm to a in 
the two different background environments shown in Fig. [2]on 
which identical internal-wave-induced sound speed perturba- 
tion fields were superimposed. In the Fig. |4]plots, dependence 
on Iq is replaced by dependence on the more familiar variable 
<p ; for an axial source this constitutes a simple stretching of 
the horizontal axis. Ensembles of AT for rays with the same 
launch angle were generated using the same technique that 
was used to generate Fig. 

Figure|4]shows that in both environments almost all of the 
structure seen in {{AT) 2 ) 1 1 2 can be attributed to the stability 
parameter a. This is because in the environments considered 
Eq. ill 2b is generally a good approximation to Eq. dlOi . Also, 
the dependence of Y,i=i on axial ray angle is somewhat 

weaker than the dependence of a on axial ray angle. We con- 
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FIG. 5: Background sound speed profiles used to construct the curves 
in Fig. 6. 

elude from Fig. |4]that travel time spreads along the timefront 
are largely controlled by the background sound speed struc- 
ture via a, and that, provided a (I) does not have too much 
structure, the dependence of Ar on a in Eq. dlOt can be taken 
outside the integral. 



B. Spreading of rays with fixed turning history and final 
depth 

We now consider a different measure of unconstrained 
travel time spread. To illustrate the generality of the results 
presented, we consider rays in upward refracting environ- 
ments (Fig. [3} in the presence of rough surface scattering. In 
such environments the action-angle form of the ray equations 
is unchanged. The action / is defined as in ijfji with the upper 
turning depth z+ = for all rays. It follows from Eqs. 0, 
neglecting d (G — W) / dr, that the travel time for one ray cy- 
cle is T(I) = 2%[I -H(I)/to(I)]. At each reflection from the 
rough surface the ray action gets modified. After N surface 
reflections and N + 1 ray cycles the travel time perturbation to 
the ray is 



AT N 



[0)(/ )P ^' 



(14) 



where Iq is the initial action of the ray, and the first of Eqs. (|3^) 
has been used. Note that Eq. (II 4i constrains the ray geometry 
and final ray depth but not the final range of the scattered rays. 
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FIG. 6: Upper panels: absolute value of the difference between per- 
turbed (rough surface) and unperturbed (flat surface) ray travel time 
after 20 surface reflections as a function of initial ray angle at the 
surface for sound speed profiles PA (a) and PB (b) shown in Fig. 5. 
Lower panels: stability parameter as a function of ray angle in each 
of the two environments. 

Figure[6]shows plots of AT^o, defined in (114-1 . vs. launch an- 
gle <p for a source at the surface [so, again, <PoCo) is a simple 
stretching] in the environments shown in Fig. [5] The rough 
surface was a single frozen realization of a surface gravity 
wavefield with a £~ 7 / 2 surface elevation wavenumber spec- 
trum with 0.02 rad m" 1 < k < 0.16 rad m" 1 , Ak = 10~ 3 rad 
m -1 , and an r.m.s. slope of 4 x 10~ 3 . Ray reflections from 
this surface were specular, but with a linearized boundary con- 
dition; the surface elevation was neglected, but the nonzero 
slope was not approximated. The same environments were 
used in Ref. |7J 

Consistent with Eq. (1141 . Fig. [6] shows that travel time 
perturbations are largely controlled by the background sound 
speed structure via a, rather than details of the rough surface. 



sion that exploits the assumed smallness of the sound speed 
perturbation 8c. The method by which the constraints are 
incorporated into the travel time perturbation estimates pre- 
sented below is different than the method that was used in 
Ref. LLP; that approach can be shown to give the same result 
that is presented below when terms of 0((AI) 2 ), neglecting 
endpoint corrections, are retained. 

The perturbation expansions presented below make use of a 
scattering model in which / following a scattered ray is piece- 
wise constant; the apex approximation is then treated as a 
special case. Away form the scattering events Eqs. are 
valid piecewise, with / making a jump at each scattering event. 
Thus Eqs. Q can be used to compute, in a piecewise fashion, 
the contribution of each ray segment to the travel time and 
range of the scattered ray. To impose the eigenray constraint, 
the total range of each scattered ray must be constrained to 
be equal to the range of the unperturbed ray. In addition, the 
turning point history, starting depth, and ending depth of the 
scattered ray must be the same as those of the unperturbed 
ray. This is accomplished by constraining the total change in 
# following the scattered ray to be equal to that of the unper- 
turbed ray. For some choices of the source-receiver geom- 
etry this procedure exactly enforces the eigenray constraint. 
In general, however, small endpoint corrections must be ap- 
plied. For simplicity, these small endpoint corrections, which 
are important only at very short range, will be neglected. 

We assume that between source and receiver the environ- 
ment can be approximated as n piecewice range-independent 
segments. Let /,■ denote the action of a scattered ray in the 2th 
segment whose horizontal extent is Ar, = At? ,•/(»(/,•), and let 
/ denote the action of the unperturbed ray. The eigenray con- 
straint requires the perturbed ray total range, r, to be equal to 
the unperturbed ray total range. This condition reads 

£ (o(it) mil) h 



IV. CONSTRAINED TIME SPREAD 

In this section we focus on deep ocean conditions and 
consider the scattering-induced broadening of an individual 
branch of the timefront. Note that this broadening is much 
smaller than the scattering-induced spreading along the time- 
front that was considered in section QUIA. We refer to the 
broadening of an individual branch of the timefront as a con- 
strained time spread because to calculate this broadening two 
constraints must be incorporated into the calculation. First, 
the rays contributing to the spread must have the same fixed 
endpoints in the (z,r)-plane. Second, the rays contributing 
to the spread must have the same turning point history, i.e. 
the same ray inclination (positive or negative) at the source, 
and the same number of turning points M (upper and lower) 
between source and receiver. Note, however, that these con- 
straints do not fix the values of the ray angle at either the 
source or receiver. 

The approach taken here to compute appropriately con- 
strained travel time spreads is based on a perturbation expan- 



which, for sufficiently small A/, = /, — 7, can be approximated 
by 

^s^-^LQs^*- <16) 

The r.h.s. term of d!6i is almost always negligible compared 
to the l.h.s. term in deep ocean environments, even at ranges 
comparable to basin scales. Then, because A#,- = o)(/,)Ar,' k 
0)(/)Ar,, the eigenray constraint reduces approximately to a 
statement that the range-weighted average of the perturbed ac- 
tion is equal to the unperturbed action, i.e. 

/W-^/iAr,-. (17) 

i=t 

This equation constrains the action history of scattered eigen- 
rays. 

The difference between the perturbed and unperturbed ray 
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travel times is 



AT 



(0(1.) 



At}; 



Ml 

(0(1) 
2 dl 2 



I A*, 



[a»(/)] 2 2dP I=r \a>)Q 

where the eigenray constraint (I16> has been used. 
Virovlyansky 14 recently derived an expression for AT 
consisting of Eq. i ll 8b [written as an integral over r using 
Ad-j = co(Ii)Ari] plus a correction term. We will discuss 
his result in more detail below. It should also be noted that 
the same symbol AT is being used to denote a constrained 
(eigenray) travel time perturbation that was used in the 
previous section to denote an unconstrained travel time 
perturbation. We believe that is obvious in all cases which 
is the relevant quantity. This choice was made to keep the 
notation simple. Also, for the same reason, we make no 
notational distinction between a theoretical estimate of Ar 
and a numerically computed Ar. 

The apex approximation is a special case of the analysis 
leading to Eqs. dl5>-(!18>. If the contributions to Ar from the 
incomplete ray cycles at the beginning and end of the ray path 
are neglected, Eq. (II 81 applies with A#; = 2n and n = N, the 
number of complete ray cycles. (For large N the neglected 
incomplete ray cycle contributions to Ar constitute small cor- 
rections to the sum of the retained contributions.) 

In Ref. [lfj a relaxed form of the apex approximation 
was considered. The transition width of the jump Air = 
had a width Ad-j, which is assumed here to be 
small, A#t < 2n. In the transition region a particular form 
h(I — #A/t/A#t) of the Hamiltonian was assumed. This 
finite width transition region was shown not to contribute 
to a range perturbation, but gives a travel time perturbation 
^A/tAtJt for a single scattering event. After N scattering 
events, assuming the transition width A#t of each is the same, 
this gives a contribution \ ~ 7o)A#t to Ar i ll 8t . while the 
eigenray constraint (I16> is unaltered. With this additional term 
( tTsl . with A&i = 2k and n = N, is replaced by 



AT = -(I N -I Q ) Adr ■ 



(19) 



For convenience we shall refer to the first and second terms on 
the r.h.s. of Eq. dl9> as Ari and A^, respectively. For steep 
rays in typical deep ocean environments Ad-j/(2n) w 0.05 
while o'A/t / CO « 0.01 . (The latter estimate is variable owing 
to variations in co'.) Thus one might expect that Ar dominates 
Ar2. For small N this is indeed the case. But under typical 
deep ocean conditions consecutive apex scattering events are 
independent (because the ray cycle distance 2%jco ~ 50 km 
exceeds the horizontal correlation length of internal waves, 
about 10 km), so A^ and Ar2 grow approximately like N l l 2 
and N, respectively. Thus Ar2 is expected to dominate A^ at 
long range. Note that when Ar? dominates Ar, this equation 
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FIG. 7: (a) Ray travel time and depth of rays with positive launch an- 
gles and which have N = 18 upper turning points at r = 1 Mm in the 
CO environment shown in Fig. 2. The solid curve correponds to rays 
in the absence of internal-wave-induced scattering. Dots correspond 
to rays scattered by internal waves, (b) Same as the left panel except 
that r = 3 Mm and N = 54. 




FIG. 8: (a) Constrained travel time spread (AT) vs. the first term on 
the r.h.s. of Eq. 1 1 91 (ATi) at r = 500 km for rays with axial launch 
angles 11° < <p < 13.5° in the CO environment with an internal- 
wave-induced sound speed perturbation superimposed, (b) As in (a) 
but vs. the second term on the r.h.s. of Eq. i 191 (ATi) and at r = 3 
Mm. 



predicts a scattering-induced travel time bias in the direction 
of sgna)'. 

These observations are consistent with the simulations 
shown in Fig. There, scattered and unperturbed (8c = 0) 
ray travel times are shown in the vicinity of two branches of 
the timefront at r = 1 Mm (N = 18) and r = 3 Mm (N = 54). 
At both ranges the source was on the sound channel axis in 
the CO environment shown in Fig. |5J launch angles for the 
rays shown are near 10° at both ranges. In that environment 
sgn co' = — 1 . It is seen in this figure that there is no indica- 
tion of a negative travel time bias at 1 Mm, while there is a 
clear negative bias at 3 Mm, consistent with Eq. d!9l >. At 
ranges sufficiently short that A^ dominates, travel time bi- 
ases of either sign can be expected as — Iq can be of either 
sign. Because In — Iq is a continuous function of launch angle, 
such biases should have a nonzero correlation scale along the 
timefront. These features are readily evident in our simula- 
tions in the presence of realistic internal-wave-induced sound 
speed perturbations. 

Quantitative tests of the correctness of Eq. J19I are shown 
in Fig. [8] There numerically computed travel time differences 
Ar are compared separately to A^ and Ar2 under conditions 
in which one of the two terms dominates the other. The con- 
strained travel time difference Ar was computed using a sin- 
gle realization of an internal-wave-induced sound speed per- 
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turbation as the difference between the perturbed ray with the 
same turning history and the same final depth as the unper- 
turbed ray. In Fig. [8]attention is restricted to rays that are suf- 
ficiently steep that the apex approximation is approximately 
valid, but not so steep that rays reflect off the surface. 

In Fig. [F^ Ar is compared to AT\ at r = 500 km; for the 
rays used to construct this plot the mean value of AT% is —0.3 
ms, which is seen to represent a small correction, on aver- 
age, to A^i. The value of Ad-j used to construct this plot is 
0.25 w 2n/25. The domain and gross distribution of AT points 
plotted in Fig. [8^ is seen to coincide with the domain and 
gross distribution of A7\ points plotted. Thus AT\ is a fairly 
good statistical descriptor of AT. The point-by-point compar- 
ison between Ar and A7J is not good, however. This is seen 
by noting that the points plotted do not cluster along the diag- 
onal line with slope unity. Thus, although AT\ has the correct 
qualitative features and appears to be a useful statistical de- 
scriptor of Ar, it is evidently a poor deterministic predictor of 
Ar. This shortcoming, we believe is attributable to the overly 
idealized form of the Hamiltonian in the apex transition region 
that was used to derive A7i . It is interesting to note that A7"i is 
the only term presented in this paper describing a travel time 
spread that is independent of and this term is the one that 
gives the poorest agreement with simulations. 

In Fig. |8j) Ar is compared to Ar? at r = 3 Mm; for the 
rays used to construct this plot the r.m.s. value of Ar is 2 ms, 
which is seen to represent a small correction, on average, to 
A^. The agreement between Ar and Ar2 in this plot is seen 
to be very good, indicating that for the rays used to construct 
this plot Ar2 is a very good predictor of Ar. Overall, for steep 
rays we have found good qualitative agreement at short range 
between simulations of Ar and A^ , and a good quantitative 
agreement at long range between Ar and Ar2 . 

Figure |9] shows plots of Ar vs. launch angle and a vs. 
launch angle at r = 3 Mm in each of the two background sound 
speed profiles shown in Fig. [2] Again, the constrained travel 
time spread Ar was computed using a single realization of 
an internal-wave-induced sound speed perturbation as the dif- 
ference between the perturbed ray travel time and the travel 
time of the unperturbed ray with the same turning history and 
the same final depth as the perturbed ray. The small gaps in 
the plot correspond to perturbed rays whose final depth lies 
outside the bounds of the portion of the timefront which has 
the same turning history as the perturbed ray. In Fig. [9]ray 
angles are not limited to the band for which the apex approxi- 
mation is expected to be valid. Thus Eq. \\9\ is not expected 
to be valid for the entire band of angles. Equation Jl 8i should 
be approximately valid across the entire band, however, so 
the trend in tt(<p ) should be approximately reproduced in the 
AT(<p ) points plotted. This is seen to be the case; variations 
in Ar are caused by variations in £" =1 (A/;) 2 A#,. The prob- 
able cause of the positive values of Ar for near-axial rays in 
the CO environment seen in Fig. (|9) will be discussed below. 
Even after fitting a smooth curve through the fluctuations seen 
in the Fig. |9]Ar plot for the C 14 environment, the peak in | Ar | 
is seen to be less pronounced than the peak in a \ . This is ex- 
pected because the scattering process leads to a local average 
over a band of adjacent a-values. With these minor caveats, 
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FIG. 9: Travel time difference Ar = 7iw — Tq at r = 3 Mm vs. ray 
launch angle <p () for an axial source, and stability parameter a vs. 
<p , in each of the environments CO (a) and C14 (b) shown in Fig. 2. 
Here, Tjw is the travel time of a ray in the presence of an internal- 
wave-induced sound speed perturbation field, and Tq is the travel of 
a ray in the same background environment, but without the sound 
speed perturbation field superimposed, whose turning point history 
and final depth are the same as those of the 7]w ray. 

Fig. |5]shows that constrained travel time spreads at long range 
are largely controlled by the background sound speed struc- 
ture via the stability parameter a. The same conclusion can 
be drawn from Fig. |8j). 

The probable cause of the positive computed values of Ar 
for near-axial rays seen in Fig. |9]in the CO environment is not 
the failure of Eq. dl9l : rather it is the failure of our ray label- 
ing scheme for near axial rays. We have implicitly assumed 
that the ray identifier ±M uniquely identifies each timefront 
branch. This assumption fails for near-axial rays. (The as- 
sumption also fails for steep rays in the vicinity of the cusps 
where adjacent timefront branches join, and whenever ©' has 
isolated zeros. We have intentionally avoided the latter situa- 
tion.) The correct way to identify timefront branches is by the 
Maslov index 8 jl which, for waves in two space dimensions, 
advances by one unit each time a caustic is touched. For near- 
axial rays the portion of the timefront corresponding to a con- 
stant value of the ray identifier consists of two adjacent partial 
branches joined at a cusp. For these rays a ray-identifier-based 
definition of Ar can lead to either undefined or nonuniquely 
defined Ar values. 

We noted earlier the quantitative correctness of Ar2 and the 
qualitative correctness of A7i . The weakest link in our analy- 
sis is clearly A^ . Virovlyansky 14 derived an alternate correc- 
tion to Ar2 whose validity, like dl 81 i. is not linked to the apex 
approximation. In our notation, this correction reads 

at? = j Q ^Ai^)—m,m^) 

- [ r ^8h(p(^,z^U). (20) 
Jo 

Unfortunately, the lack of explicit knowledge on the func- 
tional dependence of 8h on (/, r) makes the numerical eval- 
uation of the first term on r.h.s. of \20\ quite difficult — if not 
impossible. Virovlyansky argued that that term can be ne- 
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glected. We have evaluated numerically the second term on 
r.h.s. of i20\ ; this resulted in significantly poorer agreement 
with our simulations that was found using our A7i. (Note, 
however, that the parameter A#t in our expression for AT\ 
was adjusted to approximately match the simulations.) We 
expect, but have not confirmed, that the evaluation of both 
terms on the r.h.s. of i20\ would result in good agreement 
with simulations. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated three measures of travel 
time spreads for sound propagation in environments consist- 
ing of a range independent background on which a range- 
dependent perturbation is superimposed: (i) unconstrained 
spread of ray travel time along the timefront; (ii) uncon- 
strained spread of rays whose turning history and final depth 
are fixed but whose final range is not; and (iii) scattering- 
induced broadening of an individual branch of the timefront 
at a fixed location. All three measures of time spreads 
were shown to be largely controlled by a property co' of the 
background sound speed profile. Surprisingly, this is the 
same property that controls ray spreading and, hence, ray 
amplitudes 7 . We now present two arguments that provide 
some insight into why travel time spreads should be controlled 
by the same property of the background sound speed profile 
that controls ray spreading. 

The variational equations that describe how small per- 
turbations (8l,8&,8T) evolve in the extended phase space 
(/, T) are given by 




co' + 8hjj 
Ico'+l8h u 



_ Shjv _ 
I8h _/,} — 8h $ 




(21) 

where a short-hand notation for partial differentiation has been 
introduced. In general Eqs. ( 12 li and constitute a set of six 
coupled differential equations. For the special case 8h = 0, 1 
and co' are constant following trajectories and (12 H . like @, 
have a simple analytical solution. For the class of problems 
treated in this paper the nonzero sound speed perturbation 
terms in the matrix on the r.h.s. of J21t are generally much 
smaller than the co' and Ico' terms, corresponding to contri- 
butions from the background sound speed profile. Thus one 
expects that generically the dominant cause of the growth of 
(81,8$, 8T) is the background sound speed structure via co', 
rather than the small sound speed perturbation terms. Loosely 
speaking, the perturbation terms provide a seed for the growth 
of (8l,8&,8T), while subsequently growth of these quanti- 
ties is largely controlled by co'. 

Additional insight into the role played by co' is obtained 
by making a fluid mechanical analogy. Equations Q define 
a flow in the three-dimensional space (p,z,T) with velocity 
components d(p,z, T)/6r. (Recall that r plays the role of the 
independent or time-like variable in the one-way ray equa- 
tions.) Alternatively, the flow in this three-dimensional space 
can be described using action-angle variables (I,"&,T). The 



coordinates (/, #,r) behave qualitatively like cylindrical co- 
ordinates (p,9,Q, say. If we neglect the contributions from 
the sound speed perturbation term 8h, then the velocity com- 
ponents of the background flow in these coordinates are ui = 
dl/dr = 0,iif> = {d$ /dr)I = col, and u T = dT /dr = Ico H. 
The strain rate tensor, 



(22) 



where V,- denotes covariant derivative, describes how small 
elements of fluid are deformed by the flow 12 . With the 
identification (/, #, T) <-> (p, 0, £) and with the velocity field 
(ui, u$,ut) defined above, the strain rate tensor is 



\lco' \lco' 
\lco' 
\lco' 








(23) 



Although d23l > is only qualitatively correct [because of the 
qualitative connection between (/, and (p,9,Q], the 
conclusion to be drawn from this tensor is extremely impor- 
tant: deformation of small elements of the extended phase 
space (I,d-,T) by the background sound speed structure is 
caused entirely by shear (off-diagonal elements of S) and is 
quantified by the product Ico' . This behavior is consistent with 
the discussion above concerning Eq. (12 li and the growth of 
small perturbations in the extended phase space (I,t9-,T). 

Arguments similar to those given above relating to Eqs. 
(12 li and d23i were given in Ref. |7( in that study, however, 
attention was confined to ray spreading in (p,z) or, equiva- 
lently, (/,#). That problem is described by the first two of 
Eqs. Q or Q, the upper 2x2 system in Eqs. ( 12 1 i and (1231 . 
etc. In that study it was shown that ray spreading is largely 
controlled by co' . The surprising result of the present study 
is that travel time spreads are also largely controlled by co'. 
This is evident from the unconstrained and constrained travel 
time spread estimates, (1 1 Oi . ( I14> . ( I18> . and ( U9L and the more 
heuristic arguments associated with Eqs. J2 1 b and d23i . 

Inasmuch as the AET experimental observations^ pro- 
vided much of the motivation for the present work, it is note- 
worthy that the results presented here are consistent with the 
ray-based analysis of those observations presented in Ref. 
Two points relating to the results in Ref. Il0| deserve further 
comment. First, it was argued in Ref. 10 that, for moderately 
steep rays, constrained travel time spreads could be approx- 
imately estimated using AT\, the first term on the r.h.s. of 
Eq. d 1 9i . This is moderately surprising in that the AET range 
was 3.25 Mm, which is large enough that one would expect 
that the second term on r.h.s. of dl9> should be important. 
The neglected term, however, is proportional to a, which is 
unusually small (cf. Ref. 10) in that environment for the rel- 
evant band of launch angles. Thus the estimated travel time 
spread reported in Ref. [T(J is close to what one obtains using 
both terms on the r.h.s. of Eq. ( I19> . Second, it was noted 
in Ref. 10, without explanation, that in the AET environment 
constrained travel time spreads are much larger for near axial 
rays than for steeper rays. This behavior is consistent with Eq. 
( I18> and the observation 10 that in the AET environment \a\ is 
very large for the near-axial rays. 
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The results p resent ed in this paper, coupled with those pre- 
sented in Refs. I7ll4l represent an important step toward the 
development of a theory of wave propagation in random inho- 
mogeneous media (WPRIM). These studies have shown that 
both ray amplitude statistics and ray travel time statistics are 
largely controlled by the background sound speed profile, via 
ft)'. It follows that finite frequency wavefield intensity statis- 
tics should also be largely controlled by the background sound 
speed profile via (£>' . This is different from the classical treat- 
ment of the problem of wave propagation in random media 
(WPRM) which assumes that the background sound speed 
structure is homogeneous — for which ft)' = 0. Ideally one 
would like to develop a uniformly valid theory of WPRIM 
that reduces to known WPRM results in the limit ft)' — » 0. 
A more modest goal is to develop an approximate theory of 
WPRIM, appropriate for long-range underwater sound propa- 
gation, that treats the case where the wavefield intensity statis- 
tics are largely controlled by ft)'. To develop such a theory 
the results presented here and in Refs. I7ll4l have to be com- 
bined and extended. Necessary extensions are the inclusion 
of finite frequency effects (interference and diffraction) and a 
more accurate treatment of the link between sound speed per- 



turbations and perturbations to / and These topics will be 
explored in future work. 
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